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Abstract 

Open Quantum Random Walks, as developed in [I], are the exact 
quantum generalization of Markov chains on finite graphs or on nets. 
These random walks are typically quantum in their behavior, step by 
step, but they seem to show up a rather classical asymptotic behav- 
ior, as opposed to the quantum random walks usually considered in 
Quantum Information Theory (such as the well-known Hadamard ran- 
dom walk). Typically, in the case of Open Quantum Random Walks 
on nets, their distribution seems to always converges to a Gaussian 
distribution or a mixture of Gaussian distributions. In the case of 
nearest neighbors homogeneous Open Quantum Random Walks on 
Z rf we prove such a Central Limit Theorem, in the case where only 
one Gaussian distribution appears in the limit. Through the quantum 
trajectory point of view on quantum master equations, we transform 
the problem into studying a certain functional of a Markov chain on Z rf 
times the Banach space of quantum states. The main difficulty is that 
we know nothing about the invariant measures of this Markov chain, 
even their existence. Surprisingly enough, we are able to produce a 
Central Limit Theorem with explicit drift and explicit covariance ma- 
trix. In a second step we are able to extend our Central Limit The- 
orem to the case of several asymptotic Gaussians, in the case where 
the operator coefficients of the quantum walk are block-diagonal in a 
common basis. 

*Work supported by ANR project "HAM-MARK", N° ANR-09-BLAN-0098-01 
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1 Introduction 



Quantum Random Walks, such as the Hadamard quantum random walk, 
are nowadays a very active subject of investigations, with applications in 
Quantum Information Theory in particular (see [3] for a survey). These 
quantum random walks are particular discrete-time quantum dynamics on 
a state space of the form H <g> C z . The space C z stands for a state space 
labelled by a net Z d , while the space H stands for the degrees of freedom 
given on each point of the net. The quantum evolution concerns pure states 
of the system which are of the form 

After one step of the dynamics, this state is transformed into another pure 
state, 

I*') = £ Wi) ® 10 ■ 

i&Z d 

Each of these two states gives rise to a probability distribution on Z d , the 
one we would obtain by measuring the position on C zd : 

Prob({0) = |N| 2 . 

The point is that the probability distribution associated to cannot be 
deduced from the distribution associated to |\&) by "classical rules", that is, 
there is no classical probabilistic model (such as a Markov transition kernel, 
or else) which gives the distribution of in terms of the one of |^). One 
needs to know the whole state \^f) in order to compute the distribution of 

These quantum random walks, have been successful for they give rise to 
strange behaviors of the probability distribution as time goes to infinity. In 
particular one can prove that they satisfy a rather surprising Central Limit 
Theorem whose speed is n, instead of \fn as usually, and the limit distribution 
is not Gaussian, but more like functions of the form (see [5]) 

VT^a 2 (1 - \x) 
x h-» y ——^= , 

7r (1 — x 2 ) y a 2 — x 2 

where a and A are constants. 

In the article pQ is introduced a new family of quantum random walks, 
called Open Quantum Random Walks. These random walks deal with density 
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matrices instead of pure states, that is, on a state space they consider 

density matrices of the form 

To this density matrix is attached a probability distribution, associated to 
the values one would obtain by measuring the position: 

Prob({z}) = Tr ( Pi ) . 

After one step of the dynamics, the density matrix evolves to another state 
of the same form 

P Pi® , 

with the associated new distribution. 

In [1] it is proved that these Open Quantum Random Walks are a non- 
commutative extension of all the classical Markov chains, that is, they con- 
tain all the classical Markov chains as particular cases, but they also de- 
scribe quantum behaviors which cannot be described as classical stochastic 
processes. 

Though, as shown on simulations in the same article, it seems that Open 
Quantum Random Walks of infinite nets such as Z d exhibit a rather classical 
behavior in the limit, that is, their limit distribution seems to always con- 
verge to a Gaussian distribution, or to a mixture of Gaussian distributions 
(including the case of Dirac masses as particular cases of Gaussian distri- 
butions). While the quantum random walk, step by step, seems to be very 
quantum, that is, the distribution at time n + 1 has nothing much to do 
with the distribution at time n (at least it cannot be deduced from it by a 
classical process), it appears that asymptotically the quantum random walks 
becomes more and more classical. 

The aim of this article is to prove, under some conditions, a Central 
Limit Theorem for these Open Quantum Random Walks and to compute 
explicitly the characteristics of the associated Gaussian distribution: drift 
and covariance matrix. 

This article is structured as follows. In Section [2] we recall a certain 
number of notations and concepts which are very common in the context 
of Quantum Mechanics: states, density matrices, completely positive maps, 
etc. Section [3] is then devoted to presenting the general mathematical struc- 
ture of Open Quantum Random Walks and their probability distributions. 
We end up this section with a series of examples and numerical simulations 
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which illustrate our definitions and which will be covered later on by our 
Central Limit Theorems. In Section [4] we explain the Quantum Trajectory 
approach to Quantum Master Equations. This approach, which is nowadays 
very important in the study of Open Quantum Systems, gives a way for 
Open Quantum Random Walks to be simulated by means of a particular 
Banach space-valued classical Markov process. In the same section we recall 
an important ergodic property of quantum trajectories, as proved in [I]. 

The last sections are the ones where the main theorems are proved. First 
of all the main Central Limit Theorem is proved in the context of a single 
asymptotic Gaussian distribution. The proof is based on proving a Central 
Limit Theorem for a particular martingale associated to the quantum tra- 
jectories. This martingale is obtained by the usual method of solving the 
Poisson equation, which surprisingly can be implemented explicitly in our 
context, even though we do not have any information on the existence of 
an invariant measure for the Markov chain associated to quantum trajec- 
tories. Furthermore the parameters of the limit Gaussian distribution are 
explicitly obtained. We then extend the Central Limit Theorem to a context 
with several asymptotic Gaussians, but with block-diagonal coefficients for 
the Open Quantum Random Walk. We prove that, in this case, the Open 
Quantum Random Walk behaves like a mixture of Open Quantum Random 
Walks with single Gaussian, that is, up to conditioning the trajectories at the 
beginning, we get a behavior of an OQRW with a single asymptotic Gaus- 
sian. We compute several examples which illustrate the different situations 
of our theorems, we compute the associated asymptotic parameters. 

2 General Notations 

We recall here some useful notations and terminologies that shall be used in 
this article. 

All our Hilbert spaces are on the complex field and are separable (if not 
finite dimensional). For all Hilbert space H we denote by <B("H), the Banach 
space of bounded operators on "H equipped with the usual operator-norm 
that we denote by IHI^. We denote by Ci(7i) the Banach space of trace- 
class operators on "H, equipped with the trace-norm ||-|| r 

We shall be typically working on a tensor product 7-L <g> /C of (complex, 
separable) Hilbert spaces. For any <fi e /C we put {(/))%, to be the operator 

\<f>) K : H — >• U®K 
if) i — > if) ® (f> . 
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Its adjoint is the operator jc(0| defined by 

K {<f>\ : H®K, — >• ft 

x ® y i — ?• (0 , x) y . 

As a consequence of these definitions, the operator \4>)k k.{4>\ is the orthogonal 
projector onto % <8> C . 

In the case where only one Hilbert space, e.g. /C, is concerned we denote 
these operator by |0) and (0|, simply. As a consequence the operator |0)(0| 
is just the orthogonal projector onto C . 

Note that if H is a bounded operator on T-L <S> /C and if G /C, then the 
operator x;(0| ]0)jc is a bounded operator on "H. 

Recall that a density matrix p on some Hilbert space "H is a trace-class, 
positive operator such that Tr (p) = 1 . The convex set of all density matrices 
on V, will be denoted by £(H). The extreme points of this convex set are 
the pure states, that is, the rank one orthogonal projectors: 

P= 10X01, 

with G "H, || 0|| = 1. The set of pure states on T-L will be denoted by S(1-L). 

Let J\f stand for a finite or a countable set of indices. If {Ai ; i G Af} is 
a family of bounded operators on % such that 

where the convergence above is understood for the weak topology, then the 
mapping 

is well-defined, for the series is ||- (^-convergent, and the mapping preserves 
the property of being a density matrix. It is a so-called completely positive 
map. 

Note that such a completely positive map admits an adjoint map Ai* 
acting on the bounded operators on "H. More precisely, the mapping 

M*{X)=Y. A * XA ^ 

ieAf 

is a strongly convergent series and satisfies 

Tr(M(p) X) = Ti (pM*(X)) 
for all density matrix p and all bounded operator X. 
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3 Open Quantum Random Walks 



3.1 General Setup 

Let us explain here the setup in which we shall be working. It consists in 
special cases of Open Quantum Random Walks as described in p], namely, 
the case of nearest neighbors, stationary quantum random walks on Z d . Our 
presentation here is slightly different of the one of [I], for we have adapted 
our notations to the simpler context that we are studying here. 

On Z d we consider the canonical basis {ei, . . . , e^} and we put 

for all j = l,...,d. For each i G Z* d we denote by N(i) the set of its 2d 
nearest neighbors, that is N(i) = {i + ej ; j = 1, . . . , 2d}. 

We consider the space K = C z , that is, any separable Hilbert space with 
an orthonormal basis indexed by 7L d . We fix an orthonormal basis of /C which 
we shall denote by {\i)) i& d- Let H be a separable Hilbert space, it stands for 
the space of degrees of freedom given at each point of Z d . In the rest of the 
article we always assume that 7-L is finite dimensional. Consider the space 

We are given a family {A±, . . . , A 2 d} of bounded operators on H which 
satisfies 

2d 

i=i 

The idea is that the operator Aj stands for the effect of passing from any 
point i e Z rf to its neighbor i + ej. The constraint above has to be understood 
as follows: "the sum of all the effects leaving the site % is / " . It is the same 
idea as the one for transition matrices associated to Markov chains: "the 
sum of the probabilities leaving a site % is 1" . 

To the family {Ai, . . . , A 2 d} is then associated a completely positive map 
on TL, namely: 

2d 

To the family {A\, . . . , A 2 d\ is also associated a completely positive map on 
H ® fC as follows. We put 

L{ = Aj + ej}{i\ 
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for all i G Z d , all j = 1, . . . ,2d. The operator L\ emphasizes the idea that 
while one is passing from site \i) to its neighbor \i + ej) in /C, the effect on 
% is the operator Aj. It is easy to check that 

2d 

EE^ = '. 

i&L d 3=1 

where the above series is strongly convergent. Hence, there is a natural 
completely positive map on Ti <8> /C associated to these 's, by putting 

2d 

i&Z d 3=1 

for all density matrix p on H <S> K,. Recall that the series above is convergent 
in trace-norm. 

In the following, we shall be interested in iterations A4 n of Ai applied 
to density matrices of H Cg) /C. We shall especially be interested in density 
matrices on 7i <E> K. with the particular form 

p=^2p l ®\i)(i\, (1) 

where each pi is not exactly a density matrix on %: it is a positive and 
trace-class operator but its trace is not 1. Indeed the condition that p is a 
state aims to 

E Tr = 1 • ( 2 ) 

There are two reasons for such a specialization. First of all, it is easy to check 
that an application of Ad to any density matrix p on % <8> K leads to a state 
of the form Q. Secondly, those states are the states we are really interested 
in; they express no mixture between the sites, they are states which respect 
the spatial structure underlying the definition of K. 

If p is a state on H <8> K, of the form 

then a measurement of the "position" in /C, that is, a measurement along the 
orthonormal basis (|z))$ e y, would give the value |z) with probability 

p(i) = Tr (p^ . 
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After applying the completely positive map M, the state of the system H®K, 
can be easily checked to be 

MaO = £ (f>^-e^*j ®l<><i|- (3) 
iez d \j=i / 

Hence a measurement of the position in /C would give that each site i is 
occupied with probability 

2d 

p'(i) = ^Tr {A jPi _ ej A*) . (4) 

And so on, by repeatedly applying Ai to the initial state, we obtain a se- 
quence of probability measures on Z d which, in general, cannot be described 
in terms of a classical random walk. Indeed, the probability distribution at 
step n + 1 cannot be deduced from the probability distribution at step n, we 
need to know the whole states pf 1 ^ and not only their traces Tr (pf*). 

3.2 Examples 

Let us illustrate the setup above, with some examples. 

In the case d — 1, we describe a quantum random walk on Z with the 
help of only two bounded operators B and C on 7i, satisfying 

B*B + C*C = I . 

The operator B stands for the jumps to the left (it corresponds to the op- 
erator A 2 with the notations of previous subsection) and C stands for the 
jumps to the right (it corresponds to the operator Ai). 

Starting with an initial state p^ = p <g> |0)(0|, after one step we have the 
state 

p (1) = Bp B* ® | 1 > < 1 1 + Cp C* ® |1)(1| . 

The probability of presence in |-1) is Tr(BpoB*) and the probability of 
presence in |1) is Tr(CpoC*)- 

After the second step, the state of the system is 

p^ = B 2 p B 2 * ® |-2)(-2| + C 2 p C 2 * ® |2)(2| + 
+ (CB Po B*C* + BC Po C*B*) ® |0)(0| . 

The associated probabilities for the presence in |-2), |0), |2) are then 

Tr {B 2 p B 2 *), Tr {CBp B*C* + BCp C*B*) and Tr (C 2 p C 2 *) , 
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respectively. 

One can iterate the above procedure and generate our open quantum 
random walk on Z. 

As further example, take 



B 



i 



and 



C 



1 



i o 
-l i 



The operators B and C do satisfy B*B + C*C = I. Let us consider the 
associated open quantum random walk on Z. Starting with the state 



P 



(o) 



1 




®|0>(0|, 

we find the following probabilities for the 4 first steps: 



-4) 


1-3) 


1-2) 


1-1) 
1 


10) 

1 


1 + 1) 
2 


1+2) 


1+3) 


1+4) 






i 


3 


3 


3 


5 








i 


9 


5 


9 


11 


9 


10 




i 


27 


10 


27 


27 


27 


26 


27 


17 


81 




81 




81 




81 




81 



n = 
n — 1 
n = 2 
n = 3 
n = 4 

The distribution obviously starts asymmetric, uncentered and rather wild. 
The interesting point is that, while keeping its quantum behavior time after 
time, simulations show up clearly a tendency to converge to a normal centered 
distribution. Figure 1 below shows the distribution obtained at times n = 4, 
n = 8 and n = 20. 




Figure 1: An O.Q.R. W. on Z which gives rise to a centered Gaussian at the 
limit, while starting clearly uncentered (at times n = 4, n = 8, n = 20) 



A much more trivial example on Z is obtained by taking 

'1 



B 



y/p 





and C 



V^P 
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It is easy to compute the associated quantum trajectories and to show that 
they have the behavior of a random walk which goes straight to the right, 
with only one possible random jump to the left. This example will illustrate 
our Central Limit Theorem for the particular case where the Gaussian is 
degenerate. 

It is easy to produce Open Quantum Random Walks on 1? by specifying 
4 matrices N, W, S, E on H which satisfy 

N*N + W*W + S*S + E*E — I . (5) 

Then, we ask the random walk to jump from any site to the four nearest 
neighbors, following N, W, S or E, respectively. 

One can for example combine two 1-dimensional Open Quantum Random 
Walks by asking them to act on the different coordinate axis. For example, 
take 

together with 

W = v^A) (° and E = ^T^Y) (j °) , 

with a 2 + (3 2 + 7 2 = 1 and for some A e [0, 1]. 

One can obtain behaviors with a single Gaussian, as in Figure 2, with 

A = 3/4, a = 1/4, /3 = 1/4. 

The aim of the theorems to come now are to prove such Central Limit 
Theorems and to identify the elements of the limiting gaussian distribution. 

4 Quantum Trajectories 
4.1 Simulation of O.Q.R.W. 

Open Quantum Random Walks have the very nice property to admit a quan- 
tum trajectory approach, that is, a classical process simulating the evolution 
of the density matrix. This approach to Open Quantum Random Walks is 
the one that allows us to prove a central limit theorem. Let us explain here 
this approach. 

Starting from any initial state p on H <8> K, we apply the mapping M. and 
then a measurement of the position in JC, following the axioms of Quantum 
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Figure 2: An O.Q.R.W. on 1? which exhibits a single Gaussian asymptoti- 
cally ( at time n = 50) 

Mechanics. We end up with a random result for the measurement and a 
reduction of the wave-packet gives rise to a random state on H <8> fC of the 
form 

Pi <8> ■ 

We then apply the procedure again: an action of the mapping M. and a 
measurement of the position in fC. The following result is proved in pp. 

Theorem 4.1 By repeatedly applying the completely positive map M. and a 
measurement of the position on K., one obtains a sequence of random states 
on H <g) /C. This sequence is an homogenous Markov chain with law being 
described as follows. If the state of the chain at time n is = p <g) 
then at time n + 1 it jumps to one of the values 

u (n+1) = -—^AjpA/^li + eAii + eA, j = l,...,2d, 
PU) 

with probability 

P(j) = Tr (Aj pA/) . 

This Markov chain (u/™)) is a simulation of the master equation driven by 
Ai, that is, 

E [u {n+1) | co {n) ] =M(co {n) ). 

Furthermore, if the initial state is a pure state, then the quantum trajectory 
stays valued in pure states and the Markov chain is described as follows. If 
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the state of the chain at time n is the pure state \<p) ® \i) , then at time n+ 1 
it jumps to one of the values 

-~=A j \ip)®\i + e j ) , teV, 

with probability 

p{ 3 ) = \\A ] \ v )\\ 2 . 

In a more usual probabilistic language, this means that we have a Markov 
chain (p ra ,X n ) ng N with values in £{T-i) X 7L d which is described as follows: 
from any position (p, X) one can only jump to one of the 2d different values 

' Aj p Aj* , X + ej 



P(j) 



with probability 

P ti) = Tr(A jP A j 



What Theorem 4.1 says is that the law of the random variable X n coincides 
with the distribution on Z d of our open quantum random walk at time n, 
when starting with the initial state po <8> |Xo)(Xq|. 



Theorem 4.1 also says that if the initial condition is in S{Ji) £§> Z then 
the Markov chain always stays in S(H) ® Z d . 

4.2 Ergodic Property 

We now recall an ergodic theorem for quantum trajectories, as proved in [4], 
that we adapt to our context and notations. Recall the completely positive 
map on H associated to the operators A%, . . . , A 2 d ■ 

2d 

i=i 

Theorem 4.2 If (p n ,X n ) is the Markov chain obtained by the quantum tra- 



jectory procedure as in Theorem 4-1 then the sequence 

1 n 



n 

i=l 



converges almost surely to a random variable 6^ which is valued in the set 
of invariant states for C . 

In particular, if C admits a unique invariant state p^, then the above 
Cesaro mean converges almost surely to p^. 
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5 The Central Limit Theorem 



5.1 The main Theorem 

In this section we make the following hypothesis on L : 

(HI) : C admits a unique invariant state p^. 
We start with some notations. We put 

2d 

m = ^Ti(A iPoo A*) ei eR d . 
i=i 

In the following we shall denote by x ■ y the usual scalar product on Mr. We 
denote by m^, i — 1, . . . d, the coordinates of m in l d , that is = m ■ for 
i — 1, . . . , d. 

Lemma 5.1 For every I G M d , the equation 

2d 

(L - C*(L)) = A i A i (e< -l)-{m-l)I (6) 

i=i 

admits a solution. The difference between any two solutions of ^ is a 
multiple of the identity. 

Proof By definition of m we have, for every I G M, d 

2d 



/ J Tr (Ai Poo A*) ei-l = m-l 



hence 

Tr ( Poo [Y,A*A i {e i -l)-{m-l)l\ ) =0 



2d 



-2d 



We have proved that Yli=i A*Ai (ei-l) — (m-l) I belongs to {poo} • But {poo} 
is equal to Ker (I — C) 1 -, by Hypothesis (HI). Furthermore Ker (J — C) x is 
equal to the range of / — C* . We have proved that Ylf=i A*Ai (e^ - l) — (m-l) I 
belongs to the range of / — C*. This gives the announced existence. 
If L' is any other solution of (|6|) then, putting H = L — L' gives 



H-C*(H) = 0. 
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This is to say that H is an eigenvector of C* for the eigenvalue 1. By the 
hypothesis (HI), the eigenspace of £ for the eigenvalue 1 is of dimension 1. 
Hence the eigenspace of C* for the same eigenvalue is also 1-dimensional. As 
we have C*(I) = I, this means that all eigenvectors of C* for the eigenvalue 
1 are multiple of the identity. Hence if is a multiple of the identity. □ 

In the following we shall denote by L\ a solution of ^ associated to 
/ G IR d . In the case where I = e^, for some i — 1, . . . , d, we denote Li by Lj 
simply. In terms of the coordinates (Zj) of I, note that we have 

d 

Li = k Li . 

i=i 

We can now formulate our main Central Limit Theorem. 

Theorem 5.2 Consider the stationary open quantum random walk on 7L d 
associated to the operators {A\, . . . , A^d} ■ We assume that the completely 
positive map 

2d 

C(p) = J2^ P A* 

i=l 

admits a unique invariant state p^. Let (p n , X n ) n > be the quantum trajec- 
tory process associated to this open quantum random walk, then 

X n — nm 

converges in law to the Gaussian distribution A/"(0, C) in IR d ; with covariance 
matrix 

= Sij (Tr (Ai poo A*) + Tr (A+dPoo A* i+d )) - 171^+ 
+ (Tr ^ Poo A* Lj) + Tr (A, Poo A* L t ) 

— Tr (Ai +( i poo A* +d Lj) — Tr (Aj +d Poo A* +d Lj)) 

- (m, Tr (poo Lj) + rrij Tr (p^ Li)) . 

Proof Consider the Markov chain (p n ,X n ) n£ N, with values in £(1-1) x Z d , 
associated to the quantum trajectories of M.. We put AX n = X n — X n _ 1; 
for all fief and we consider the stochastic process (p n , AX„) ngN . which is 
also a Markov chain, but with values in £("H) x {ex, . . . , e2d}- Its transition 
probabilities are given by 

[ otherwise, 
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for all i, j G {1, . . . , 2d}. 

We are given a fixed I G M d and we wish to write a Central Limit Theorem 
for (X n ■ Z) n£ N- Our first step is to find a solution to the so-called Poisson 
equation, that is, we wish to find a function / on £(%) x {ei, . . . , t2d} such 
that 

(I-P)f(p,x)=x-l-m-l. (7) 



is 



Lemma 5.3 Equation |?p admits a solution which 

f(p,x)=Tr(pL l )+x-l. (8) 



Proof [of Lemma 5.3 



If we define / by 

f(p,x) = Ti(pLi)+x-l, 

we get 

(2d 
J2 r ^(ApA*L l )+ 
i=i 

2d 

+ J2Tr(A iP A*)ei-l 

i=i 

= Tr (p (^Lt - C\U)) - ^A*A iei -l^ +x-l 
= —m ■ I + x ■ I . 

That is, the function / is a solution of the Poisson equation. D[of Lemma] 

The second step of the proof consists in carrying the problem of our 
central limit theorem to a central limit theorem for a martingale. 
With the help of the Poisson equation, we have 

n 

X n -l-n{m-l)=X -l + J2 (( X k ~ *k-i) -m)-l 

k=i 

n 

= X -/ + ^(/-P)/(p fc ,AX fc ) 
k=i 

n 

= X -l + J2 (f(Pk, AX fc ) - Pf( Pk -i, AX fc _0) 



k=2 

+ f( Pl ,AX 1 )-Pf( Pn ,AX n ) 
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We put 

n 

M n = Y, f(Pk, ^X k ) - P/(p fc _!, AI H ) . 

k=2 

Clearly (M n ) n > 2 is a centered martingale, with respect to the filtration (J r n )„>2, 
where T n = a{(p k ,X k ) ; k < n}. We put 

R n = X ■ I + f(p u AXx) - Pf(p n , AX n ) . 

We claim that (|-R n |)neN* is bounded. Indeed, by Equations ^ and Q we 
have 

Pf(p n ,AX n ) = Ti(p n L l )+m-l 
and \Ti(p n Li)\ is bounded independently of n by 

llPnlli Halloo = ll^lloo • 

This means that the term R n has no contribution to our central limit theo- 
rem. It is thus sufficient to obtain a central limit theorem for the martingale 
(M n ) ng N*. We recall the form of the Central Limit Theorem for martingales 
that we shall use here. 

Theorem 5.4 (cf [2J, Theorem 3.2 and Corollary 3.1) Let (M n ) n6N be 
a centered, square integrable, real martingale for the filtration (J r „) ne N- If, 
for all e > 0, we have the following convergences in probability: 

1 n 

lim - E [(&M k ) 2 l ]A M k \>eVE I -Ffc-i] = (9) 

k=l 

and 

1 n 

lim TE [(AM,) 2 | JT fc _J = a 2 (10) 

k=l 

for some a > 0, then M n j ^Jn converges in distribution to a A^(0, a 2 ) distri- 
bution. 

As a third step of our proof we shall prove that (M n ) ra > 2 satisfies the 
property (|9]). We have 

AM k = f(p k , AX k ) - P/(p fc _x, AI H ) 

= Tr (p k Li) + AX k - l-m-l-Tn (p fe _i ■ L t ) . 
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In particular AM k is bounded independently of k for 

|AMfe| < \\pk\\i Halloo + \\^ X k\\ ¥\\ + \\ m \\ ¥\\ + Ibfc-llll Halloo 

< 2 ||L,|L + 11211 + IMI IUII . 



The condition (|9| is then obviously satisfied as ^\AM k \>e^n vanishes for n 
large enough. 

The fourth step of the proof consists in computing the quantity 

E[(AM fc ) 2 | F k _ x ] , 



in order to verify that Condition ( 10 ) is satisfied. We have 

AM k = Tr (p k L x ) - Tr {p k _ x L x ) + {AX k - m) ■ I 

so that 

(AM,) 2 = Tr (p h Uf - Tr {p h _ x U) 2 

- 2 Tr {p k _ x U) [Tr (p k L x ) - Tr { Pk _ x L x ) + (AX k - m) ■ I) 
+ (AX k ■ I - m ■ I) 2 + 2 Tr (p k L t ) (AX k -l — m-l). 

We denote by T±, T 2 and T 3 , respectively, the three lines appearing in the 
right hand side above. The term E[Ti | J^k-i] is equal to 

E[Tr (p k L t ) 2 | F k „ x \ - Tr (p k L x ) 2 + Tr {p k L{f - Tr (p k ^ L x ) 2 . 

The term E[Tr (p k L x ) 2 \ J^k-x] ~ Tr (p k L x ) 2 is the increment of a martingale 
(Y n ) and it is bounded independently of k (using the same kind of estimates 
as for \R n \ above). Hence Y n /n converges almost surely to 0. 

The term Tr (p k L x ) 2 —Ti (p k -i L x ) 2 , when summed up to n gives Tr (p n L x ) 2 — 
Tr (pi L x ) 2 and hence converges to when divided by n. 

The term E[T 2 | J^-i] clearly vanishes for it makes appearing the condi- 
tional expectation of the increment of the martingale (M n ). 

Note that here appears a key point in our proof: all the quadratic terms 
in p k disappear in the limit; this is crucial for otherwise it would have been 
impossible to handle them without information on the invariant measure of 
the Markov chain (p n ). 
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We finally compute E[T 3 | Tk-i)- We get 

E[T 3 | J- fc _i] = E [(AX, ■ I) 2 - 2(m ■ /)(AX fe • I) + (m ■ lf+ 

+2 Tr ( Pk L{) (AX k ■ I - m ■ I) \ T k _ x ] 

2d 

= Tr (A p fc _! A?) [(e, • /) 2 - 2(m • I) (e, • /)] + 

i=l 

2d 

+ 2^Tr (Ai^A*^) (e 4 -/-m-/) + (m-0 2 

= Tr A *i A * {e l -l-m-lf + 

+ 2A* L l A i (e i - 1 - m ■ 1) 
We put 

2d 

Ti = J2 A * A i (ei-l-m-l) 2 + 2A* L { A, (a ■ I - m ■ I) . 

i=i 

Putting everything together, by the fact that Y n /n converges to and by 



the Ergodic Theorem |4.2[ we get that 

1 



^E[(AM fe ) 2 |.7V a ] 

'" k=3 

converges almost surely to 

of = Tr (poo Ti) . 

The fifth and last step of the proof consists in rewriting the variance of 
in order to make the covariance matrix C appearing. We have 

2d 2d 



Ti = J2 A *M e * ■l-m-l) 2 + 2j2A*L l A i (e i -l-m-l) 

i=l i=l 
2d 2d 

= J2 A *M e i ■ If - 2(m • /) A *M e i -l) + (m- 1) 2 + 

i=l i=l 
2d 

+ 2 J2 A* U A{ ei ■ l) - 2(m • l)£*(Lt) . 



i=i 
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Hence, this gives 

2d 

Tr ( Poo r,) = Tr ( A i P°° A *)( e * ■ l ? ~ 2 ( m • l f + (rn ■ lf+ 

2d 

+ 2 Tr (A Poo ^ Lj)(ei • - 2(m • Z) Tr (£( Poo ) L,) 



i=l 

2d 2d 



= -(m • Z) 2 + ^ Tr (A Poo A*)(e, • /) 2 + 2 ^ Tr (A Poo A* L,)(e i • Z) 
i=i i=i 
-2(m-Z)Tr( Poo L ; ). 

This gives 

d d 

^ = " J2 m * m M + Yl l i ( Tr P°° A i) + Tr Poo A *i+d)) + 

i,j=l i=l 
d 

+ 2 ^ Z i Z j (Tr (A Poo A* Lj) - Tr ( A l+d Poo L,)) 

d 

- 2 ^ ZiZj-mj Tr ( Poo %) • 
This proves that 

d 

°f = ^ Cij , 

where the matrix a is the one given in the theorem statement. The theorem 
is proved. □ 

5.2 The one dimensional case 

The one dimensional case is a useful one, we make simpler in this case the 
formulas we have obtained above. 

In the case where the dimension is d — 1, there are only two jump oper- 
ators A\ and A 2 , which satisfy 

A\A ± + A* 2 A 2 = I . 

We have 



m = Tr (A, Poo A\) - Tr (A 2 Poo A* 2 ) . 
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In dimension 1 there is only one operator Lj, the operator Li, which we 
denote here by L simply and which is solution of 

L - C*{L) = A\A X - A* 2 A 2 -ml = 2A\A 1 - (1 + m)I . 

Finally, following the theorem above, we have 

a 2 = Tr (A lPoD Al + A 2Poo A*) - m 2 + 

+ 2 Tr \{A lPoo A\ - A 2Poo A* 2 ) L] - 2m Tr ( Poo L) 
= 1 - m 2 - 2m Tr ( PoD L) + 2 Tr [{A lPoa A\ - A 2Poo A* 2 ) L] 
= 1 - m 2 - 2m Tr ( Poo L) + 2 Tr ( PoD L) - 4 Tr [(A 2Poo A;) L] 
= 1 - m 2 + 2(1 - m) Tr ( Poo L) - 4 Tr [ Poo A* 2 LA 2 ] 
= 1 - m 2 + 4 (Tr (A 2Poo ^) Tr ( Poo L) - Tr ( Poo A;LA 2 )) , 

or else 

a 2 = l-m 2 + A(Ti (p^AlLA,) - Tr {A lPoa A\) Tr ( Poo L)) . 



5.3 Examples 

We shall now explore several examples in order to illustrate our Central Limit 
Theorem. Let us first start with two examples on Z. The example 

B = ^ (o alld = ^(-i \ 

that we mentioned earlier falls in the scope of our theorem for it admits a 
unique invariant state 

Poo = 

In particular we have 

m = Tr (C Poo C*) - Tr (B Poo B*) = . 

We recover here that the limit Gaussian distribution is centered, as was 
observed in the simulations above. 

is 



The operator L, given by Lemma 5.1 



1/5 -1 
3 1-1 



+ XI ■ 



This gives 

a 2 
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Let us compute the case of our trivial example on Z obtained by taking 



In that case the unique invariant state is 

Poo = 



1 


We find m = 1 in that case, which is compatible with the behavior we 
described for this example. 

The operator L in this case is 



-2 




+ XI. 



This gives a 2 = 0. We recover that the asymptotic behavior of this open 
quantum random walk is degenerate, with drift +1. 

Let us end up this illustration with the 2-dimensional example mentioned 
in Subsection EX 



1 A A 1 (1 0\ TTr 1 A) 1\ „ 1 A 

We find a unique invariant state 



"=2 1 ' S =2 -1 1 ■ W =S 1 ■ £ 410 



_ 17 
Po ° - 33 I 16 

The average is 

29 -1\ 

m = I — — , — — 



132' 132 J 

The two solutions of Equation ^ are then 

68(16+^) \ / 8(l6+-v/i4) 







_ 3993 r _ 3993 

j(l64Vl4) 8(756+17^) I ' I 8( 1 6+v / Ti) 4(-57+4v / li) 



3993 3993 / \ 3993 3993 

and we find the following covariance matrix 

'0.675 0.008^ 



( ' 0.008 0.211 



approximately . 
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6 The Block-Diagonal Case 



6.1 The Main Theorem 

The Central Limit Theorem proved above does not concern the case where 
C admits several invariant states. This is typically the case when the asymp- 
totic behavior shows up several Gaussian contributions. The proof we have 
obtained above does not adapt to the general case. However, there is one 
situation, with several Gaussians which we are able to treat. Let us describe 
it now. 

Consider the operators A±, . . . , A 2 d satisfying 

2d 

Ea*a = /, 

i=i 

as previously. We now assume that there exists a decomposition 

U = E 1 © E 2 © . . . © E N 

of "H into orthogonal subspaces such that all the A^s are block-diagonal with 
respect to this decomposition. That is, 

Ai(Ej) C Ej 

for alH = 1, . . . , 2d, all j — 1, . . . , N. This hypothesis is denoted by (HI') in 
the rest of this section. We denote by Pj the orthogonal projector onto Ej. 
Note that the condition above is equivalent to 

PjAi = A{Pj 

for alii = 1, . . . , 2d, all j = 1, . . . , N. We put 

A^P = AiPj = PjAiPj . 

In the same way we denote by the completely positive map associated 
to the operators (A^)f^ v On each subspace Ej we have 

2d 2d 
i=l i=l 

If p is a density matrix on % we put 

P U) = I'll'') ■ 
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Let p be a density matrix and P p the law of the Markov chain (p n , X n ) n > 
obtained as previously, by the quantum trajectories associated to the matri- 
ces Ai, starting with the initial state p. Recall that 

/ A tPn A* \ 

with probability Tr (Aip n A*). 
We put pi ] =Tr(P jPn ). 

Lemma 6.1 The process (pn^) n >o ^ s a martingale for the filtration 

F n = o-((p k ,X k ) , k<n) . 

Proof We have 

[p% | Fn] = Tr ( P 3 A iPn A *) 

i 

i 
i 

= ^Tt (A*A iP ^) 

i 

= Tr(p(f))=^'). 

□ 

As is a martingale we can consider the associated Girsanov trans- 

form (that is, the /i-process). We define Pp ^ to be the law on the trajectories 
which is given, on the length n trajectories by 

Po 

where P n is the law on the trajectories with length n. In other words 



p U) p 



where p$ = Mmpffl. 
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Proposition 6.2 Under the law the sequence 

JJ) 

Pn v 



Tr (p^) 



n>0 



has the law of the quantum trajectories associated to the family of operators 
(-<4j^)i=i,...,2d an d starting from the state 

Proof The sequence p$ = Tr (P jPn ), n e N, is a function of (p„). The chain 
(p n ,X„) under P( J ) is thus a /i-process of the initial chain for the harmonic 
function p^\p) = Tr(Pjp). We thus have that (p n ,X n ) is a Markov chain 
under with transition probabilities: 

' „ Aip„A* 

Pn+1 



with probability 



Tr (A iPn A*^ 
X n+ \ = X n + Ci 



^Tr (A iPn A*) 

Pn 



But we have 



% Tr(A iPn A*) = ^f^, Tr (A iPn A*) 

Tr(P jPn ) 

Tr(4 j) pj j) 4 jr ) 
Tr(P jPn ) 

We see that the transition probabilities only depend on the component p\P. 
If we consider the sequence 

Ji) 

■Jj) _ P n 
Pn Tr(p^) 



we have 

X n+ i = X n + Ci 



24 



with probability 

Tr (A^A?*) 
Tr(#) 

This exactly means that the sequence X n ) n > under has the law of 
the quantum trajectories associated to the family (A^)^. □ 

We now make the following hypothesis. 
(H2) Each of the mappings admits a unique invariant state 
We then put = (mf \ . . . , m^j) where = Tr (A k p^A* k ). 
(H3) The m^'s are all different. 

Under these hypotheses we have the following result. 

Theorem 6.3 Under the hypotheses (HI 7 ), (H2) and (H3) we have the fol- 
lowing properties. 

1) For allj = 1,...,N, 

F p [UmpW = 1] = p { j) = 1 - P p [limp^ = 0] , 

that is, the vector p n = {pf , . . . ,Pn) converges to (0, ... , 0, lj, 0, . . . , 0) with 
probability p^ (note that YljPo^ = ^-)- 

2) Conditionally to limpn = 1 (that is, under the measure 

P p [-| limp^ = 1] = Wf) 

we have that (]M\X n ) has the law of the quantum trajectories associated to 
the family of matrices {A\ )f=i- In particular, under this conditional law, 
the process 

(X - nm®) 



n 



converges in distribution to the Gaussian distribution J\f (0, C^V where 



is given by the same formula as in Theorem 5.2 but for the family (A^). 



Note that the theorem above concretely means that the quantum trajec- 
tories in that CclSG £11*6 du mixture of Open Quantum Random Walks of the 



form of Theorem 5.2 The associated stochastic process can be obtained as 
follows: with probability p^ the process (X n ) follows the law of the Open 

Quantum Random Walks with associated matrices A± and then satisfies the 
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corresponding Central Limit Theorem with mean and covariance matrix 

Proof We know that under Pj^ the sequence (f$\x n ) has the law of 
the quantum trajectories associated to the family (A^). As the mapping 
£W admits a unique invariant state we also know that if we consider N n {i) 
to be the number of jumps e$ made by the quantum trajectory up to time n, 
then we have 



lim —N n (i) = m 
n 

00 



CO 



almost surely for the measure ¥] 

We also know that (p$) is a martingale; it is furthermore non-negative 
and bounded, hence it converges almost surely and in L 1 to a limite 
This means that the support Q^i) f P^ is given by 

n«) = {pg > 0} . 

If the ra^'s are all different then the measures Pj^ 
a consequence the sets are all different and finally p, 
probability 1. 

Note that we have 



p are all singular. As 



(?) 

'oo 



or 1 with 



P 



for 



F p [.| limj#> = 1] 
P„[.| limrf> = l] = 



P. 



limp 



0") 



but P 



Jo) _ i 

Poo — J- 



Po^ for (pn') is a martingale. 



The conclusion now is a direct consequence of the Central Limit Theorem 
established for the chain (X n ) but now associated to the family (A^)^ v 

□ 
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